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Abstract 

The quantum mechanical motion of a relativistic particle in a non- 
continuous spacetime is investigated. The spacetime model is a dense, 
rationale subset of two-dimensional Minkowski spacetime. Solutions 
of the Dirac equation are calculated using a generalized version of 
Feynman's checkerboard model. They turn out to be closely related 
to the continuum propagator. 

Keywords: Dirac equation, Feynman checkerboard, discrete space- 
time 

PACS: 03.65.P, 03.30, 04.20.G 

1 Introduction 

Continuity and a metric of Lorentz signature are fundamental properties 
associated with spacetime in the framework of Special and General Rela- 
tivity. Though intuitive, the assumption of continuity is debatable because 
of quantum mechanical arguments. Discrete models of spacetime have been 
investigated thus by various authors ||, ||, |, | . 

The modern idea that space could be discrete goes back to Riemann |^] 
who reflected about a natural measure of space. More recent attempts to 
consider spacetime as discrete either aim at resolving divergence problems in 
Quantum Field Theory at a fundamental level or seek to reconcile General 
Relativity and Quantum Theory to form a unified theory. 

In this paper we investigate possible effects, which the hypothetical 
discreteness of spacetime could have on Quantum Theory. To this end. 
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Feynman's path integral approach to Quantum Tehory ?] provides an 
ideal framework since it naturally accounts for spacetime properties. For 
continuous spacetime, the path integral formalism has been shown to be 
completely equivalent to ordinary wave mechanics. However, for a non- 
continuous spacetime one should not necessarily expect equivalent findings. 
To gain insight into the possible differences between quantum mechanics 
in the framework of common continuous spacetime and a non-continuous 
spacetime model, we seek solutions for the Dirac equation. 

As has been pointed out by Feynman and Hibbs Q, the retarded prop- 
agator of the 1 + 1 dimensional Dirac equation 

id-^/dt = -ia^d-^/dx - a^-^ (1) 

(with units c = Ti/m = 1) can be obtained from a random walk model 
in which the particle motion is restricted to movements either forward or 
backward at the speed of light. We consider this model in the framework 
of a particular non-continuous spacetime model and demonstrate by an 
explicit calculation that the model yields common solutions of the Dirac 
equation. 

The paper is organized as follows: Section 2 gives a short review of 
relevant properties of the non-continuous Minkowskain spacetime model de- 
scribed in In section 3 we calculate solutions of the Dirac equation. 
Section 4 summarizes and discusses the results. 



2 Non-Continuous Minkowskian Spacetime 

The two dimensional spacetime model considered here is defined by the 
subset M of R2 

M = {t,x] = {^{p^ + q^),^{p^-e)]<lB?, (2) 
m m 

n,m,p,q e Z \ {0} 

In the following, n,m,p,q may take any value of the indicated range and 
are not to be considered as fixed. Elements of M correspond to spacetime 
points with temporal and spatial coordinates {t,x). M is rational, hence 
countable, i.e. of cardinality while R'^ is of cardinality Hi. With respect 
to this difference we denote M as non-continuous. 
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The set M is invariant under the transformations 



'):M^M, (3) 
4>{s) = 




2pq \ -{p^-q^) P^ + q^ 
Z\{0}. 



The parameter 



2 2 

X p — q 



V = v(t, x) = - = ^ (4) 

^ ' t p^ + q^ ^ ' 

corresponds to the velocity defined as ratio of space x and time t displace- 
ments measured from the origin of the coordinate system. As illustrated 
by the second line of (^), the transformation (f) maps rational spacetime 
coordinates (t, x) onto rational ones {t',x'). 

It is straightforward to verify that (i) the elements f £ (p niap onto M, 
that (ii) that (p with respect to matrix multiplications o possesses group 
structure and that (iii) {(p, o) is a subgroup of the common 1 + 1 dimen- 
sional Lorentz group (using natural units). According to these properties 
and the non-continuity of M we denote M as non- continuous Minkowskian 
spacetime. In P] it has been shown that in the framework of the d-space 
formalism |l^, 11 1 the set M corresponds to a spacetime model which is 
nowhere diffeomorphic to but possesses the key properties of the com- 
mon macrophysical spacetime. In the same context M may be considered 
as a massless solution of the Einstein equations. Note that in contrast to 
other discrete spacetime models (e.g. [Q, |5|, and references therein) which 
intrinsically violate relativistic covariance and introduce an observationally 
not supported minimal length, the spacetime model considered here exhibits 
a generalized form of covariance. Since the set M is dense in R^, it doesn't 
introduce a minimal length as well. 

The following sections are based on the spacetime model We will 
consider light-cone coordinates (r, I) = {^{t+x), ^{t — x)) in order to simplify 
calculations. In light-cone coordinates M becomes 

M* = {r,l} = {-p^-q'}cK' (5) 
m m 

m,n,p,q £ Z \ {0}. 
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To end this section we note that since the set M is dense, it is possible 
to formulate differential equations without referring to a calculus of finite 
differences. 



3 Quadratic Checkerboard Model 

In 1^ Feynman and Hibbs described a model for the quantum mechani- 
cal motion of a relativistic electron. Solutions of the 1+1 dimensional Dirac 
equation are obtained by summing over all possible particle trajectories con- 
sisting of movements either forward or backward at the speed of light. As- 
suming natural units c = fi/m = 1, the motion of the particle corresponds 
to a sequence of straight path segments of slope ±45° in the x-t plane. The 
retarded propagator tp^^{x,t) of the Dirac equation is obtained from the 
limiting process (e.g. [@, [l3|| ) 

i^s^{x,t) = lim As^ie) Ns^{R){ief (6) 

R>0 

N is the number of segments with constant length e = t/N between the start 
point (which is assumed to be the origin of the coordinate system) and the 
end point {x,t) of the path. R denotes the number of bends while Ns^{R) 
stands for the total number of paths consisting of N segments with R bends. 
The indices 7 and 5 correspond to the directions forward or backward at the 
path's start and end points, respectively, and refer to the components of ip. 
As'y{€) accounts for the appropriate normalization. 

Since the early 80'ies, the checkerboard model of relativistic particle 
motion has been subject to anew interest (see e.g. [|l3|, U, |l5|, |l^, p^, p^). 
In Ijl^ it has been observed that paths with fixed start and end points have 
R — 1 degrees of freedom, i.e. the last bend of a path is fully specified by 
the location of the R — 1 preceding bends and the end point of the path. 
As well it has been shown that the normalization constant As.y{€) = 1 if 
only the R — 1 bends enter the calculation which actually define the path 
of the particle. The general intuition behind this is the idea that the origin 
of the divergence problems of Quantum Field Theory could be related to 
an over-specification of the theory. Here, the corresponding interpretation 



of the checkerboard model and the calculation scheme described in |19] are 
adapted. 

In the following we demonstrate by an explicit calculation that expression 
(0) when generalized to account for the key properties of M yields solutions 
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of the Dirac equation. Before starting with the calculation, some remarks 
about the limiting process @ and the terminology are in place. First observe 
the structure of (^). ijjs-^ix^t) is defined in the limit — > oo where 
is the number of path segments a path has. The length of each of these 
path segments is e = t/N . Bends of the particle trajectory occur only at 
boundaries of path segments. This is equivalent to consider a particle moving 
on a rectangular spacetime lattice with equal spacing length (or 'resolution') 
e. Consider now the expression 

Ns,{R){ief . (7) 

i?>0 

which is subject to the limit oo. For a given A^, it is weighted sum 

of the number Ns^{R) of possible paths with A^ segments and R>0 bends 
linking the start and the end point under consideration. The complex weight 
(amplitude) each path is contributing is (ie)^, i.e. each bend of a path 
contributes (ie). As has been demonstrated (e.g. j^, |l3|), taking the limit 
N ^ oo yields solutions of the Dirac equation. 

When adapting Feynman's checkerboard model to the discrete space- 
time M we will follow as closely as possible the scheme outlined above. Two 
differences have to be taken into account, the first of which being a gener- 
alization. First, it is supposed that path segments can be of non-uniform 
length en 7^ const, and that the contribution of each bend to the overall 
amplitude of a path is proportional to the segment length e^, immediately 
preceding the bend. This leads to the generalized form 

R 

i^5^{x,t)= lira As^{e)y2Ns^{R){ifr\en (8) 

^0 n=l 

of (^). Second, only particle paths on M will be considered. These path are 
characterized by bends occurring at spacetime points (t, x) satisfying 

t p^ + q^' ^ ' 

In terms of M* , the light-cone representation of M, this condition becomes 

I q'^ 

where (r, /) = (^(t + x), i(t — x)). Points of the spacetime lattice should ac- 
count for restriction ([To|). We achieved this by considering quadratic space- 
time lattices whose lattice points {jr,ji) S have, up to scaling, spacetime 
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lightcone coordinates (r, Z) = {jr,jf)- The spacetime coordinates of lattice 
points of this kind naturally account for condition (^). Figure |l| shows 
an example for a path with = 8 segments. With respect to the modi- 
fied properties of the spacetime lattice, the model is denoted as quadratic 
checkerboard model. 

Prom the definition of the quadratic spacetime lattice (or directly from 
Figure ||) it is evident that path segments are of non-uniform length and 
are depending on the lattice location {jr,ji)- We only will be interested in 
the length of paths segment preceding a bend and bends of a path to the right 
and to the left will be considered separately in the following. Observing that 
the length of a path segment preceding a bend to the right (left) is uniquely 
specifying by the lattice coordinate jV {ji), the indexes r, / can be dropped. 
Then, the length of a path segment immediately following a bend can be 
written as 

ej = {2j - l)eo (11) 

The constant eo accounts for the scaling of the path. It depends on the 
end point (t, x) and on the number of segments a path has. For an 
explicit expression of eo in terms of the number of segments to the left 
and right, see equation (p^. Note that eo scales not just the length of paths 
segment followed by a bend but the whole quadratic spacetime lattice. From 
this it follows immediately that for N ^ oo the spacetime lattice under 
consideration becomes a rectangular subset of M. 

In line with the generalization (|8|) of Feynman's checkerboard model sup- 
pose now that each bend defining a path on M contributes an amplitude 
proportional to the length of the path segment immediately preceding the 
bend. Again, bends to the right and left are considered separately and the 
following notation is adopted: The index n enumerates the segments of a 
path followed by a bend to the right (left). j„ stands for the lattice coor- 
dinate jr {ji) of the segment and the set indicates the path segments, 
after which bends to the right (left) occur. For a schematic example consider 
Figure The length ej,^ of the path segment followed by the n*'* bend to 
the right (left) according to (|ll|) becomes 

e,„ = {2jn - l)eo (12) 

and the corresponding contribution of the n*'' bend to the overall amplitude 
is 

^jn = iejn = (2jn - l)eo (13) 
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The total amplitude of a path contributed by bends to the right (left) is 
given by the product 

n 

Having established this notation, we start with the evaluation of (|8|). As 
in the 'linear' case (|19[|) a path with R bends that starts with a positive 
velocity (i.e. to the right) and ends with a negative velocity (i.e. to the left) 
consists of exactly (i? — 1)/2 + 1 bends to the left and {R — 1) /2 bends to the 
right. The {R — l)/2 bends to the right can occur after an arbitrary path 
segment to the left. {R — l)/2 of the (-R — l)/2 + 1 bends to the left occur 
in the same manner after path segments to the right while the additional 
bend to the left must occur after the last segment to the right. Let P denote 
the total number of path segments to the right (+) and Q those to the left 
(— ). In total, the path has N = {P + Q) segments. The contribution of the 
R+ = {R- l)/2 bends to the right is 

p-i 

V'-+(i?+)= (2ji-l)-...-(2j.-l)-(ieo)^^ (15) 

Next consider the situation where the path consists of a large number of 
segments to the right, i.e. P ^ 1 or equivalently — > oo. This limiting 
process corresponds to the spacetime lattice becoming dense. For P ^ 1, 
V' i-(-^^) is approximated by 

V-+(P+) « 75^ E (ji-...-J„) -2^^(160)^^ (16) 

(17) 

(18) 
(19) 

The contribution of the R^ = (i? — l)/2 + l bends to the left is calculated 
similarly. The additional bend (occurring after the last segment to the right) 
does not enter the calculation since a path is fully determined by the R — 1 
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bends to the right and left, respectively. We find 

2(ii--i)^ig ^(R--l) /02\ 

^-.iR-) ^ [l] (20) 



Q2(i^--l)(ig^)(fi--l) 

{R--iy. 



(21) 



To get the total contribution of the R bends to the left and the bends 
to the right, expressions ip V' i-(-R ) multiplied yielding 

^_+(i?) = (i6o)^-^T777?^SMi2 (22) 



[((^-l)/2) 



which represents the contribution of all paths with exactly R bends that start 
to the right and end to the left. Finally, the summation over all possible 
paths that start to the right and end to the left is preformed by summing 
V'_+(i2) over the total number of bends R. This leads to 

where 

^0 = (24) 

accounts for the proper scaling of the lattice. Notice that for ^ oo 
the scaling factor eo 0. Defining the classical velocity of the particle 

V = Ax/ At = x/t = (P2 _ Q2)/(p2 + g2) g^^^ pQ ^ ^p2 ^ Q2)/2-y 

where 7 = 1/ Vl — v"^. Finally h becomes 

- P)!]2 

= Mt/l) (26) 

A similar calculation of V'h leads to the same result. This can be seen 

from interchanging the roles of P, Q, R" and R^ . 

For ■)/'++) the number of bends to the right and to the left is R/2 for each 
direction where R is even. However, the path again is defined by R'^ = i?/2 
bends to the right and R" = R/2 — 1 to the left. Thus, 

= i:/o)''-' w2).W2-l)! 



even . 
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_ ,Py._,.,mf:^ (30) 

= (31) 

With P/Q = (x + t)/{t^ — x^Y^"^ the component becomes 

^^^ = -Sl±^^J,(th). (32) 
A similar calculation for tp leads to 

V^__=i^^7Ji(t/7). (33) 

which completes the calculation. 

To relate the components ij^^^ to the Dirac equation (||) consider its 
explicit representation with 




In this representation, tpi and ij)2 defined as 

V'l = f 5++ V V'2 = ( I (35) 




are two common independent solutions of the Dirac equation (|T]). Since M 
is dense in R^, the differential form of (||) does not pose a problem even 
though M is non-continuous. 

The four components tps-y are closely related to the matrix elements of 
the retarded continuum propagator of the Dirac equation (see e.g. p3[| ). 
Despite this there is no full equivalence and a physical interpretation of 
in terms of propagator components is not straightforward. In the last section 
we will come back to this and related questions in more detail. 

As a side remark note that for the applied calculation scheme it is es- 
sential to consider only those bends which actually define the path of the 
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particle. While it is straightforward to 'renormalize' Feynman's original 
model if all bends are considered, this is not the case here. Taking into ac- 
count all bends, the normalization 'constant' As^{€o) for the checkerboard 
restricted to M turns out to not only depend on eo but also on 6, 7 and the 
classical velocity v. 

4 Discussion and Conclusions 

It is counter-intuitive that Feynman's spacetime based path integral ap- 
proach to quantum theory yields in the case of the discussed example in- 
triguing similarities between the common continuous and the discrete space- 
time model M. The solutions ip obtained from the generalized checkerboard 
model and the common retarded continuum propagator of the Dirac equa- 
tion are of the same form. 

However, despite the same functional form, the components ipS"/ are not 
fully equivalent to the retarded continuum propagator. Reworking the out- 
lined calculation with physical units h/m ^ 1 shows the dimensionality of 
ips-y to be lenght^ which is different from the dimensionality lenght~^ of the 
retarded continuum propagator. With respect to the original checkerboard 
model, the difference is related to two origins: First, taking into account only 
R — 1 bends means accounting for one factor me less. Second, Feynman's 
original checkerboard model is 'renormalized' by dividing components like 
e.g. (p3|) and (|3^) by 2e in order to yield the retarded continuum propagator. 

The dimensionality lenght^ of ips^y may be seen as suggesting a picture 
consistent with the cardinality of M. Profoundly taking into account the 
cardinality Hq, a continuum- like propagator is inconsistent with M. This is 
due to the fact that because of its cardinality there exists no measure on 
M. As a consequence, the concept of integration is not available. Instead, 
one can expect (infinite) sums to play the role of integration. This suggests 
that the components of a discrete counterpart of the continuum propagator 
(i.e. ipS'y) should have dimensionality lenght^ . Fully clarifying the physical 
meaning of in the framework of the discrete spacetime model M is not 
straightforward. One is lead into number theoretical problems that are not 
easy to tackle. 

Other open questions remain. For example: To what extend does the 
path integral approach capture the relation between spacetime and quan- 
tum mechanics? This question is raised not only by the present results but 
already by Feynman's original checkerboard model. After all, according to 
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its formulation, bends occur only at spacetime points with rational coordi- 
nates. As well, the velocity spectrum of the particle is rational. In this sense, 
Feynman's checkerboard model is of cardinality and thus discrete in the 
same manner as M is. To our knowledge this is a fact that the literature 
has not paid attention to. 

In ||2^ it has been shown that the solution of the Dirac equation can be 
understood in terms of bit strings detached from any notion of spacetime. 
The present results seem to point into the same direction, namely, that there 
is (at least in the special case of the Dirac equation) only a loose connection 
between spacetime properties and quantum mechanics, except probably for 
the signature of the metric. 

Other open issues are related to the spacetime model M itself. The ratio 
x/t interpreted as classical velocity of the particle (starting from the origin) 
implies a discrete velocity spectrum 

2 2 

W = ^2^' P,9GZ\{0}. (36) 



What are the implications of the non-continuous spectrum (3£)? Progress 
in this and other directions again is complicated by the fact that M is of 
cardinality Hq. Finally, a question not addressed here is where the bends of 
a path occur, that contribute most to the overall amplitude. Because of the 
structure of M it is not evident that the locations are distributed as in the 
continuous case. 

In summary, it has been shown that Feynman's path integral approach 
restricted to the subset M of recovers common solutions of the Dirac 
equation. The particle path on M is limited to bends at positions i with 
coordinate ratios Xi/ti = {p^ — q^) / {p^ + q^) where p,q £ Z \ {0}. The 
bends of paths on M thus occur at a subset of the rational space and time 
coordinates only. It should be noted that both the characteristics of M and 
the discrete velocity spectrum (|3^) are invariant under the Lorentz subgroup 
(f). In particular, bends of paths occur only at rational space and time 
coordinates and the velocity remains rational in all coordinate systems. In 
this sense the model exhibits a generalized form of Lorentz invariance missing 
in Feynman's original checkerboard model and most other lattice models 
of spacetime. This demonstrates that there is a discrete, two-dimensional 
spacetime model of cardinality accounting for some of the key properties 
of special relativity and quantum mechanical features of relativistic particle 
motion. 
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Figure 1: A path on M between the origin of the coordinate system and a 
given end point {t,x). The path has a total oiN = P + Q = 8 segments, 
where P = 5 segments are going to the right and Q = 3 segments to the 
left. A total oi R = i?+ + R~ = 5 bends occur, i?+ = 2 to the right and 
R~ = 3 to the left. Note that the last bend to the left is fully determined by 
the end point {t, x) and the location of — 1 = 4 bends, i.e. by the i?+ = 2 
bends to the right and the first R~ — 1 = 2 bends to the left. 
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